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Abstract 

Within a hadron-string dynamical transport approach (HSD) we investigate the 
attenuation of high transverse momentum (j>±) hadrons as well as the suppression 
of 'near-side' and 'far-side' jets in Au + Au collisions at invariant energies y/s = 200 
GeV and y/s = 62.4 GeV in comparison to the data available from the Relativistic 
Heavy-Ion Collider (RHIC). From our transport studies we find that a significant 
part of the high j>j_ hadron attenuation seen experimentally can be attributed to 
inelastic interactions of 'leading' pre-hadrons with the dense hadronic environment. 
In addition, we also show results of 'near-side' and 'far-side' angular correlations of 
high p± particles from Au+Au collisions at y/s = 200 GeV and y/s = 62.4 GeV 
within this (pre-)hadronic attenuation scenario. It turns out that the 'near-side' 
correlations are unaltered - in accordance with experiment - whereas the 'far-side' 
correlations are suppressed by up to ~ 60% in central collisions. Since a much larger 
suppression is observed experimentally for these reactions in central reactions we 
conclude that there should be strong additional (and earlier) partonic interactions 
in the dense and possibly colored medium created in Au+Au collisions at RHIC. 
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1 Introduction 

The phase transition from partonic degrees of freedom (quarks and gluons) to 
interacting hadrons is a central topic of modern high-energy physics. In order 
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to understand the dynamics and relevant scales of this transition - as well 
the dynamical formation of hadrons - laboratory experiments are presently 
performed with ultra-relativistic nucleus-nucleus collisions. Estimates based 
on the Bjorken formula [1] and/or calculations within hadron-string-dynamical 
transport approaches [2,3,4] indicate that the critical energy density for the 
formation of a quark-gluon plasma (QGP) of 0.7 to 1 GeV/fm 3 [5] should by 
far be exceeded in the initial phase of central Au+Au collisions for a couple 
of fm/c at top SPS and RHIC energies [6]. However, the task still remains to 
unambiguously identify the formation and properties of this new phase. 

Whereas there are a couple of tentative signatures for an early partonic phase 
such as the high transverse pressure seen experimentally [7,8] as well as the 
large elliptic flow of hadrons at moderate and high momenta [7,9], a unique 
'smoking gun' could not be identified so far. In order to find representative 
signatures one has to look for early 'hard' probes that dominantly test the 
first few fm/c of the relativistic heavy-ion collision. Hence, the measurements 
of 'hard' jets should offer a direct access to the dynamics in the early stage 
when the (possibly deconfined) matter is very dense. In fact, the PHENIX 
[10] and STAR [11] collaborations have reported a large relative suppression 
of hadron spectra for transverse momenta p± above ~ 3—4 GeV/ c which might 
point towards the creation of a QGP, since this suppression is not observed in 
d+Au interactions at the same bombarding energy per nucleon [12,13,14] and 
thus has to be attributed to final state interactions. Furthermore, an almost 
complete suppression of 'far-side' jets has been reported for central Au+Au 
collisions [16] which is also not seen in d+Au reactions [17]. 

These observations - at first sight - appear to be compatible with the picture of 
partonic energy loss in a dense colored medium [18,19,20,21,22,23,24,25,26,27,28,29,30,31], 
however, it is not clear presently to which extent these suppression phenom- 
ena might be due to ordinary hadronic final-state-interactions (FSI) [32] since 
(in-)elastic collisions of (pre-) hadronic high momentum states with the bulk 
of hadrons in the late fireball might also contribute significantly to the atten- 
uation of the high transverse momentum hadrons at RHIC [33,34,35]. Hence, 
the hadronic attenuation has to be addressed in more detail before conclusions 
on perturbative QCD effects in a deconfined QGP phase on the materializing 
jets can be drawn. In fact, as shown in Refs. [35,36] a major fraction of the 
high p± hadron attenuation in Au+Au collisions at \fs = 200 GeV can be 
described by the (dominantly) inelastic interactions of 'leading' pre-hadrons. 

In order to specify the latter findings in Refs. [35,36] we briefly recall the 
concept of 'leading' (pre-hadronic) and 'secondary' (ordinary) hadrons in the 
string-hadron transport approach: In a high energy nucleon-nucleon collision 
two (or more) color-neutral strings are assumed to be formed; the string ends 
are defined by the space-time coordinates of the constituents which are denoted 
as 'leading' quarks. The latter are initially colored but assumed to pick up 
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almost instantly a quark (or antiqiuark) from the vacuum to build up color 
neutral 'leading' pre-hadrons [37]. These pre-hadronic states are distinguished 
from 'secondary' hadrons that arise from the further fragmentation of the 
strings and consist solely of quarks and antiquarks (diquarks and antidiquarks) 
created from the vacuum at some later time. The total time that is needed for 
the fragmentation of the strings and for the hadronization of the fragments 
is denoted as formation time and taken as r/ ~ 0.8 fm in the HSD transport 
approach. Due to time dilatation the formation time tf in any reference frame 
is then proportional to the Lorentz 7-factor, 

t f = -fr f . (1) 



We stress that the formation time is a sensible parameter and has been fixed 
in [38] by the rapidity distribution of charged particles at SPS energies in the 
range 0.6-0.8 fm/c. Though Tf is no longer a free parameter we will perform 
also calculations for the jet suppression in Section 4 for different Tf ranging 
from 0.2 to 0.5 fm/c. In principle, other scenarios for modeling the string 
fragmentation might also be investigated, e.g. formation times depending ex- 
plicitly on the hadron species or those extracted from the LUND fragmentation 
scheme (cf. Ref. [39]); however, we here constrain to the simple model intro- 
duced in Ref. [38] and delay an investigation of alternative prescriptions to 
future studies. 

In the HSD transport approach, furthermore, it is assumed that hadrons, 
whose constituent quarks and antiquarks are created from the vacuum in the 
string fragmentation, do not interact with the surrounding nuclear medium 
within their formation time tf and thus do not create early pressure. On 
the other hand, for the leading pre-hadrons, i.e. those involving quarks (an- 
tiquarks) from the struck nucleons, a reduced effective cross section ai ea d is 
adopted during the formation time tf and the full hadronic cross section later 
on. The effective cross section oiead is fixed in line with the additive constituent 
quark model [35] (cf. also Section 3). 

As shown in Ref. [35] for central Au+Au collisions at s/s = 200 GeV hadrons 
with transverse momenta larger than ~ 6 GeV/c predominantly stem from the 
string ends and therefore can, in principle, interact directly with the reduced 
cross section o~i ea d ■ Moreover, we have to stress that in the HSD approach 
the formation of secondary hadrons (e.g. in Au+Au collisions) is not only 
controlled by the formation time tj, but also by the energy density in the 
local rest frame: hadrons are not allowed to be formed if the energy density 
is larger than lGeV/fm 3 , which is about the critical energy density for a 
phase transition to the QGP in thermal equilibrium [5]. The interaction of the 
leading and energetic (pre-)hadrons with the soft hadronic and bulk matter 
is thus explicitly modeled to happen only for local energy densities below 
that cut, i.e. pre-hadronic interactions with fully formed hadrons as well as 
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interactions between formed hadrons are delayed until the energy density has 
dropped below lGeV/fm 3 in the expansion phase! Note, however, that the 
highest energy densities are achieved in the center of the 'fireball' and that 
the formation of hadrons starts at the surface where substantially lower energy 
densities are encountered. 

We recall that the interactions of leading (pre-) hadrons are also important to 
understand the attenuation of hadrons with high (longitudinal) momentum in 
ordinary 'cold' nuclear matter. The studies in Refs. [39,40] (see also [41]) have 
shown that the dominant final-state-interactions of the hadrons with maxi- 
mum momentum - as measured by the HERMES Collaboration [42] - are 
compatible with the concepts described above. These independent investiga- 
tions impose severe constraints on hadron attenuation in a medium that might 
possibly be attributed to a QGP phase (cf. also Ref. [29]). 

In this work we concentrate on the transverse momentum dynamics and es- 
pecially on the very high momentum tail of the hadron spectra. In order to 
describe these high p ± spectra, we use the Pythia v6.2 event generator [43] 
for nucleon-nucleon collisions and the HSD transport model [38,44,45] for the 
space-time localization of the events and further propagation in time including 
final-state interactions (cf. Refs. [35,36]). 

Our study is organized as follows: We start in Section 2 with an analysis 
of high pj_ hadron events from pp reactions - generated via PYTHIA v6.2 
- with respect to correlations in the azimuthal angle tp employing different 
momentum cuts for the hadrons. Section 3 gives a very brief reminder of 
the HSD approach and the actual 'input' to the transport calculations. The 
calculated results for jet quenching in Au+Au collisions at y/s = 200 GeV 
for different centrality classes are presented in Section 4. The sensitivity of 
the jet suppression to the formation time r/ will be demonstrated, too. In 
Section 5 we will show the results of calculations for Au+Au collisions at 
y/s = 62.4 GeV for different centrality classes since this system is presently 
under investigation experimentally. Apart from a comparison of the high p± 
attenuation of charged hadrons with the data from the PHOBOS collaboration 
[46] we will also present our predictions for the 'near-side' and 'far-side' angular 
correlations at this energy. Section 6 summarizes our present investigations. 



2 Jet correlations in pp reactions 

Before coming to the more subtle issue of jet quenching in Au+Au collisions 
at RHIC energies, it is instructive to have a look at the jet topology in pp reac- 
tions. As mentioned before, we employ the event generator PYTHIA v6.2 which 
has been adjusted with respect to the (kj_) parameter to properly describe the 
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transverse hadron spectra at \/s = 200 GeV (cf. Fig. 1 in Ref. [35]). As a 
trigger particle we choose a charged high p± hadron in the rapidity interval 
\y C m\ < 0.5 with momentum 4 GeV/c < p T ^ 19 < 6 GeV/c as in Ref. [15]. The 
correlated particles then are selected in the interval 2 GeV/c < p s ^ c < p 7 ^ 19 
and \y cm \ < 1.0. The resulting angular correlation 

. x 1 dN 



with respect to the azimuthal angle of the trigger particle is displayed in 
Fig. 1 (upper left part) and shows two distinct peaks which are attributed 
to 'near-side' and 'far-side' jet correlations since the maxima at ip=0 and 
ip = 7r correspond to back-to-back emission. Furthermore, the comparison to 
the data from the STAR Collaboration [15] shows that the angular correlations 
measured for pp reactions at \fs = 200 GeV are well reproduced by the event 
generator Pythia v6.2. 
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Fig. 1. The correlation in the azimuthal angle <p for pp reactions at -y/s = 200 GeV 
for different cuts on the trigger particle and associated particles from the event 
generator Pythia v6.2. The data in the upper left part are taken from Ref. [15]. 
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The structure and relative height of the two peaks changes when selecting 
momenta 8 GeV/c < p^ 19 < 10 GeV/c for the trigger particle and 4 GeV/c 
< p s ^ c < p 7 ^ 19 for the other charged hadrons (cf. Fig. 1, lower left part). 
In case of the higher momentum trigger particle the angular distribution of 
associated hadrons (from jet fragmentation) becomes narrower; the 'far-side' 
peak gets larger since the number of associated hadrons in the peaks are lower 
than for the cuts in the upper left part of Fig. 1. The asymmetry in the height 
of the peaks is due to the fact that the 'trigger particle' is not counted in the 
'near-side' angular correlation. Apart from the sensitivity to the momentum 
cuts demonstrated above the height of the peaks also depends on the rapidity 
cut for the associated particles. When gating on a larger rapidity interval 
| Ay | = 2 (r.h.s. of Fig. 1) the 'near-side' peak roughly stays the same whereas 
the 'far-side' peak increases by almost a factor of two. 

A more detailed information is provided in the upper part of Table 1 where the 
probability of finding 0, 1, • • •, 3 further hadrons in the 'near-side' and 'far-side' 
peak is given for 4 GeV/c < p T ^ ig < 6 GeV/c and 2 GeV/c < p^ c < p 7 ^ 19 and 
\y cm \ < 1-0. These numbers demonstrate, that in less than 10% of all events 
a further particle is found in the given rapidity and p± interval and that the 
conditional probability to find any further particle roughly decreases by an 
order of magnitude. These low conditional probabilities should be kept in 
mind when comparing azimuthal angular correlations from Au+Au collisions 
(cf. Sections 4 and 5). We note additionally that 'jet structures' only show up 
after averaging over a large event sample. 

It is, furthermore, instructive to have a look at the rapidity distribution of 
charged particles when gating on a trigger hadron in the rapidity interval 
|ycm | < 0.5 with momentum 4 GeV/c < p™ 9 < 6 GeV/c. The correlated 
particles in the 'far-side' peak - selected in the interval 2 GeV/c < p s ^ c < 
pTrig _ sriow a ra ther wide distribution in rapidity as seen from Fig. 

2. Here, the upper solid line corresponds to the correlated charged hadrons 
in the 'far-side' peak (color: red) whereas the lower solid line corresponds to 
the correlated charged hadrons in the 'near-side' peak (color: green) which 
is localized in rapidity around the trigger particle. The dashed lines give the 
rapidity distribution for the correlated charged hadrons in the 'far-side' (upper 
dashed line, color: red) and 'near-side' peak (upper dashed line, color: red) 
which correspond to leading particles, i.e. where the hadron carries at least 
one of the constituent quarks of the initial protons. Note that in pp collisions 
two strings are formed with in total four leading particles. For the momentum 
cuts adopted the leading particles show up dominantly at \y\ ~ 4, i.e. in the 
rapidity regime of the projectile/ tar get remnants. 

Consequently, the hadrons produced and/or formed by the fragmentation of 
jets do not occupy a narrow rapidity window close to the trigger hadron but 
are widespread in rapidity at least in the 'far-side' peak (with even a small 
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Table 1 

The conditional probability to find in addition to the trigger particle 0,1,. .. 3 parti- 
cles in the 'far-side' and 0,1,. . . 3 particles in the 'near-side' peak in pp (upper part) 
and central Au+Au collisions (lower part) at y/s = 200 GeV. The trigger conditions 
are: 4 GeV/c < p[ ri9 < 6 GeV/c and 2 GeV/c < p s ± ec < p[ rig and \y cm \ < 1.0. 

dip at midrapidity) . This wide distribution also explains why the probability 
of observing a correlated hadron together with a trigger hadron in roughly the 
same rapidity interval |Ay| m 1 and high momentum is rather low (cf. upper 
part of Table 1). 



3 Ingredients of the transport model 

We employ the HSD transport model [38,44,45] for our study of p+A, d+A 
and A+A collisions. This approach takes into account the formation and mul- 
tiple rescattering of formed hadrons as well as unformed 'leading' pre-hadrons 
and thus incorporates the dominant final state interactions. In the transport 
approach we include nucleons, A's, N*(1440), N*(1535), A, S and S* hy- 
perons, S's, S*'s and f2's as well as their antiparticles on the baryonic side 
whereas the 0~ and Y~ octet states are included in the mesonic sector. Inelas- 
tic hadron-hadron collisions with energies above ^fs ~ 2.6 GeV are described 
by the Fritiof model [47] (employing Pythia v5.5 with Jetset v7.3 for the 
production and fragmentation of jets [48]) whereas low energy hadron-hadron 
collisions are modeled in line with experimental cross sections. We mention 
that no explicit parton cascading is involved in our transport calculations 
contrary to e.g. the AMPT model [49] or explicit parton cascades [50]. 
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Fig. 2. The rapidity distribution of correlated hadrons when gating on a high p± 
hadron in the interval 4 < p T ™ g < 6 GeV/c for \y\ < 0.5 and 2 < p^ c < p T ™ 9 . 
The 'near-side' contribution is also peaked at midrapidity (lower solid line, color: 
green), however, shows only a small content of 'leading' particles (lower dashed line, 
color: green). On the other hand, the rapidity distribution of correlated hadrons in 
the 'far-side' peak (upper solid line, color: red) is widespread with maxima around 
\y\ rs 4 dominantly due to leading hadrons (upper dashed line, color: red). 

A systematic analysis of HSD results and experimental data for central nucleus- 
nucleus collisions for 2 AGeV to 21.3 ATeV has shown that the spectra for 
the 'longitudinal' rapidity distribution of protons, pions, kaons, antikaons and 
hyperons are in reasonable agreement with available data. Only the pion ra- 
pidity spectra are slightly overestimated from AGS to SPS energies [51,52,53] 
which implies, that the maximum in the K + /tt + ratio at 20 • • -30 AGeV - seen 
in central Au+Au (Pb+Pb) collisions [54] - is missed. Furthermore, there 
are more severe problems with the dynamics in the direction transverse to 
the beam. Whereas the pion transverse momentum spectra are rather well 
described from lower AGS to top RHIC energies the transverse momentum 
slopes of kaons/ antikaons are clearly underestimated above ~ 5 AGeV in cen- 
tral Au+Au collisions. In Ref. [53] this failure has been attributed to a lack of 
pressure generation in the very early phase of the heavy-ion collisions which 
also points towards a new state of very strongly interacting matter in the very 
early phase of central Au+Au (or Pb+Pb) collisions. In addition, the ellip- 
tic flow of high p± hadrons is missed substantially [35] which also suggests 
strongly interacting matter in the very early phase. 

Inspite of the deficiences pointed out the overall reproduction of the exper- 



8 



imental spectra is sufficiently realistic such that we can proceed with more 
exclusive probes that are produced and propagated in the background of the 
expanding string/hadron matter generated in relativistic nucleus-nucleus col- 
lisions. 

3. 1 Perturbative treatment of high p± hadrons and jets 

For the production and propagation of hadrons with high transverse mo- 
mentum (> 1.5GeV/c) we employ a perturbative scheme as also used in 
Refs. [2,55] for the charm and open charm degrees of freedom and in Ref. 
[35] for high p± hadrons. 

The initial conditions for the production and also subsequent propagation of 
hadrons with moderate to high transverse momentum (> 1.5 GeV/c) are incor- 
porated in the HSD approach by a superposition of pp collisions described via 
PYTHIA [43] which serve as the basic input. The partons (from PYTHIA) 
fragment to hadrons - according to the JETSET-part of PYTHIA - with 
the vacuum fragmentation functions. We mention that medium-modified frag- 
mentation schemes might have to be employed for nucleus-nucleus collisions, 
however, we will only use the vacuum fragmentation functions in this study 
in order to find out, if this minimal scheme also applies for the heavy-ion 
case. Together with the 4-momenta of all resulting 'particles' we store the 
constituent quark content and the relative weight W; L of each event. 

The weight Wi is given by the ratio of the actual production cross section 
divided by the inelastic nucleon-nucleon cross section, e.g. 

W, = ^gfef^ . (3) 

a NN W s ) 

All 'particles' with one or more constituent quarks are denoted as 'pre-hadrons' 
while those without are denoted as 'secondary' hadrons. The pre-hadrons start 
propagating in space-time from the space-time position of the individual in- 
teraction vertex x l v , which is well defined from the transport calculation, and 
begin to interact immediately with a fractional cross section according to the 
constituent quark model (see below). The 'secondaries' are inserted in the 
dynamical calculation after their formation time tf (in the actual reference 
frame) at a space-time position which is displayed from the vertex x l v by 
Vj ■ tf with Vj denoting the velocity of the 'secondary' hadron j. The veloc- 
ity Vj is fixed in the calculation by the momentum and energy of the finally 
formed hadron j, which is an approximation. In practice the major fraction 
of 'secondaries' shows up with velocities close to 1 (in units of c) such that 
alternative modelings have no major effect on the final results for attenuation. 
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We then follow the motion of the perturbative hadrons within the full back- 
ground of strings/hadrons by propagating them as free particles, i.e. neglecting 
in-medium potentials, but compute their collisional history with baryons and 
mesons or quarks and diquarks. For reactions with diquarks we use the corre- 
sponding reaction cross section with baryons multiplied by a factor of 2/3. For 
collisions with quarks (antiquarks) we adopt half of the cross section for col- 
lisions with mesons and for the leading pre-hadron (formed) baryon collision 
a factor of 1/3 is assumed. The elastic and inelastic interactions with their 
fractional cross section are modeled in the HSD approach in the same way as 
for ordinary hadrons with the same quantum numbers via the FRITIOF model 
[47] (including Pythia v5.5 with Jetset v7.3 for the production and frag- 
mentation of jets [48]). This concept - oriented along the additive quark model 
- has been proven to work rather well for nucleus-nucleus collisions from AGS 
to RHIC energies [51,52,53] as well as in hadron formation and attenuation 
in deep inelastic lepton scattering off nuclei [39,40]. We stress again that the 
latter reactions are important to understand the attenuation of pre-hadrons 
or ordinary hadrons with high momentum in 'cold' nuclear matter [37]. Our 
studies in Ref. [39,40] have demonstrated that the dominant final state interac- 
tions of the hadrons with maximum momentum, as measured by the HERMES 
collaboration [42], are compatible to the concepts described above. This also 
holds for antiproton production and attenuation in proton-nucleus collisions 
at AGS energies [56]. 



3.2 Simulation of the Cronin effect in pA, dA and A A collisions 

As known from the experimental studies of Refs. [57,58] an enhancement of 
the high transverse momentum particle cross section from proton-nucleus col- 
lisions - relative to scaled pp collisions - is already observed at SPS and ISR 
energies. This 'Cronin effect' is probably related to an increase of the average 
transverse momentum squared (k\) of the partons in the nuclear medium. 
One may speculate that such an enhancement of (k\) is due to induced initial 
semi-hard gluon radiation in the medium, which is not present in the vacuum 
due to color neutrality [25,37]. Since the microscopic mechanisms are beyond 
the scope of our present analysis, we simulate this effect in the transport ap- 
proach by increasing the average (k±) in the string fragmentation with the 
number of previous collisions Ap rev as [35,53,59] 

(kl) = (kl)Jl + aN pTCV ) . (4) 

The parameter a ~ 0.25 — 0.4 in (4) has been fixed in comparison to the ex- 
perimental data for d+Au collisions at RHIC [13,14] in Ref. [35]. As discussed 
in our previous analysis [35] the Cronin effect has an essential impact on the 
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shape of the nuclear modification factor 



Raa(p±) 



1/N%%* d 2 N AA /dydp ± 




(5) 



(N con ) /< clas d 2 a pp /dydp ± 



used experimentally to quantify the in-medium effects from d+Au or A+A 
collisions relative to pp reactions. In (5) the quantity (N co u) denotes the num- 
ber of binary nucleon-nucleon collisions - directly available from the transport 
calculation - while 0"™ elas stands for the inelastic nucleon-nucleon cross section 
(known experimentally). 



4 Results of Transport calculations for Au+Au collisions at ^Js = 



We recall that the detailed analysis of high p± hadron attenuation at the 
top RHIC energy of y^i = 200 GeV in Ref. [35] has demonstrated: i) the 
(pre-) hadronic final state interactions are able to approximately reproduce 
the high p±_ suppression effects observed in d+Au and Au+Au collisions; ii) 
the interactions of formed hadrons after a formation time tf ~ 0.8 fm/c are 
not able to explain the attenuation observed for transverse momenta p± > 6 
GeV/c, whereas the shape of the ratio Raa (5) in transverse momentum p± 
reflects the presence of final state interactions of formed hadrons in the 1-4 
GeV/c range; iii) the attenuation by (pre-) hadronic final state interactions 
is insufficient to explain the large suppression seen experimentally in central 
Au+Au collisions. Since our transport calculations are identical with those 
presented in [35] we directly continue with the new analysis on jet suppression. 

Fig. 3 shows our results for Au+Au collisions at y^s = 200 GeV for different 
centrality classes with respect to the 'near-side' and 'far-side' angular corre- 
lation of charged hadrons. We have gated on high p± particles in the interval 
(p± ig = 4 ... 6 GeV/c) and accumulated further charged hadrons in the inter- 
val p± = 2GeV...p^" 3 < 0.7). In addition to the angular correlations 
from Au+Au collisions we have added the results from pp reactions (dashed 
lines) in comparison to the data from STAR for collisions [15] (full dots) for 
each centrality bin to demonstrate both the actual statistical uncertainty due 
to different event samples for the varying centralities (dashed lines) and the 
actual attenuation generated by the final state interactions in the transport 
calculation (solid lines, color: green). 

We find that - when gating on high p± hadrons (in the vacuum) - the 'near- 
side' correlations are close to the 'near-side' correlations observed for jet frag- 
mentation in the vacuum (cf. lower part of Table 1) for all centralities. This is 
in agreement with the experimental observation [9,15]. Furthermore, for the 
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'far-side' correlations we get a ~60% reduction for the most central collisions 
but not an almost complete disappearance of the 'far-side' jet as indicated by 
the experimental data [9,15]. The latter data for central Au+Au collisions are 
displayed in the upper left part of Fig. 3 by the open squares. The suppression 
for 5-20 %, 20-45% and 45-65% centralities in the 'far-side' peak is lower than 
for the 5% most central reactions, but not dramatically different. 
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Fig. 3. 'Near-side' and 'far-side' jet-like correlations from HSD for p+p 
and central Au+Au collisions at midrapidity for p T ™ 9 = 4 ... 6 GeV /c and 
p± = 2GeV/c. . .p T ^ 19 {\y\ < 0.7) at different centralities. Data for p+p (full dots) 
are from [15]. The dashed lines show our input to the transport calculations for 
each centrality class, while the solid lines show the result including the final state 
interactions (color: green). The experimental data [9,15] for central Au+Au colli- 
sions are displayed in the upper left part by the open squares (color: red). Note, 
however, that the experimental trigger is p 7 ^ 19 = 4 . . . 6GeV/c and further charged 
hadrons have been accumulated in the interval p± = 2 GeV . . .p T ™ 9 for |Ary| < 0.5 
with respect to the trigger particle. Thus the trigger conditions differ slightly from 
those in the calculations. 



In order to become more quantitative, we show in Fig. 4 the integral over the 
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'far-side' peaks (cf. Fig. 4) 
Sl% 2 dtp C(A<p) 



AA 



LAA 



Slf dp C(Aip) 



(6) 



VP 



as a function of the number of participating nucleons A part normalized to the 

respective integral for pp reactions. The HSD calculations (full dots) indicate 

a rather fast drop with A part and almost become constant for A part > 150. 

We note that experimentally the decrease of a very similar ratio appears to 

continue more strongly with A part [9] (open squares with error bars in Fig. 4). 

1 /s 

For orientation the shaded area (color: green) shows an attenuation ~ A part , 
where the relative strength has been fixed by either the calculated result at low 
or high Ap ar t, respectively. Apparently, this power law in A part is too strong 
in comparison to the transport calculation (full dots). 




400 



A 



part 



Fig. 4. The angular integral over the 'far-side' correlation (6) for Au+Au at ^/s = 200 
GeV from HSD (full dots with statistical error bars) as a function of the number of 

1/3 

participating nucleons A part . The shaded area shows an attenuation ~ A ' t , where 
the relative strength has been fixed by the calculated result either at low or high 
Apart, respectively. The solid line results from the simple estimate (7) when plotted 
versus A par t and fixed in size at high A par t- The experimental data [9] are displayed 
in terms of the open squares (color: red) where the statistical and systematic errors 
have been added quadratically. Note, however, that experimentally the associated 
charged hadrons have been accumulated in the interval p± = 2 GeV . . -p^ 19 for 
| Ar]\ < 0.5 with respect to the trigger particle. Furthermore, the angular integrations 
have been carried out in a slightly larger angular range. 



We note that a rather good reproduction of the HSD results is obtained when 
comparing to the quadratic average of the long and short axis of the almond 
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shaped interaction zone - in x, y direction - as a function of impact parameter 



with R L = y R 2 — 6 2 /4, R s = R — 6/2 and R denoting the radius of the 
Au-target in hard-sphere geometry. The resulting suppression is displayed in 
terms of the solid line in Fig. 4 - fixed in absolute magnitude at high A part 
- using the relation A part {b). In view of the agreement between the rough 
estimate (7) and the numerical results of the HSD calculations we conclude 
that the suppression in the transport calculations is essentially determined by 
a linear geometrical length scale, i.e. by the average jet propagation length 
in the almond shaped reaction zone (cf. Ref. [60]). Note in addition that the 
relation A part {b) is model dependent for very peripheral reactions (small A part ). 

Since the suppression of 'far-side' jets from pre-hadronic final state interactions 
in our transport calculations is less than 60% at all centralities, there should 
be an additional suppression mechanism at least for mid-central and central 
collisions that is not included in the present transport approach. Since the 
dynamics at intermediate and late times are expected to be properly described 
by HSD only the very early phase of the nucleus-nucleus collision should be 
failed. 

We recall that in HSD (as well as UrQMD [3,4]) there are dominantly collisions 
of pre-hadrons (string ends) in the early collision phase which neither generate 
enough early transverse pressure [53], elliptic flow [52,61,62] or attenuation of 
high p± hadrons [35,36]. These general findings are in line with those seen in 
the suppression of 'far-side' jets in Figs. 3 and 4: the medium of strings (pre- 
hadrons) is not interacting strongly enough to explain an almost complete 
quenching of 'far-side' jets in central Au+Au collisions at y/s = 200 GeV. 
Nevertheless, about 60% (in central Au+Au collisions) is still a considerable 
amount of attenuation that should be missed in parton scenarios based on 
perturbative QCD cross sections. 

As mentioned before the hadron formation time r/ is a sensible parameter of 
the transport approach although fixed by hadron rapidity distributions [38]. 
Nevertheless, one might argue that this formation time is simply not short 
enough and an alternative 'fixing' might yield the proper jet suppression. 
A similar analysis e.g. has been performed in Ref. [62] with respect to the 
strength of the elliptic flow t> 2 of charged hadrons, where shorter formation 
times Tf lead to an increase of the interaction rate at short times and conse- 
quently to a rise of the elliptic flow. To this aim we have performed calculations 
for 5% central Au+Au collisions at y/s = 200 GeV introducing an indepen- 
dent formation time 77 for the jets and their fragmentation products without 




(7) 
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Fig. 5. The same as in Fig. 3 for most central collisions (upper left part), but now 
for formation times Tf = 0.5 fm/c (l.h.s.) and 0.2 fm/c (r.h.s.) instead of the default 
value of 0.8 fm/c in Fig. 3. 

changing the formation time of the bulk 'non-perturbative' hadrons. The re- 
sults of these calculations for the azimuthal angular correlations are shown in 
Fig. 5 for Tf = 0.5 (l.h.s.) and 0.2 fm/c (r.h.s.) by the solid lines (color: green). 
In fact, the attenuation of the 'far-side' peak increases with decreasing forma- 
tion time. However, also the 'near-side' peak now shows a reduction by the 
final state interactions as seen from a comparison of the dashed and solid lines. 
This reduction of the 'near-side' peak is still very modest for 77 = 0.5 fm/c but 
reaches a level of ~ 50% for ry = 0.2 fm/c which can clearly be ruled out by 
the present data. In short: a sizeable reduction of the formation time - though 
leading to a more substantial suppression of the 'far-side' peak as well as to 
an increase of the elliptic flow of charged hadrons [34,62] - is incompatible 
with the experimental observations for central Au+Au collisions. 



Without explicit representation we mention that a decrease in the energy den- 
sity cut for hadron formation (cf. Ref. [35]) from 1 GeV/fm 3 to 0.7 GeV/fm 3 
leads only to a small reduction of the 'far-side' jet suppression. This also holds 
for a moderate increase of the cross section for the pre-hadronic states. These 
results are easy to understand: The 'leading' pre-hadrons in central collisions 
are attenuated by a factor of about 15 (cf. Fig. 8 of [35]) with the standard 
cross sections, such that these particles are practically erased. Furthermore, 
only 70-80% of high p± hadrons are pre-hadrons (cf. Fig. 3 in [35]). We point 
out that a larger cross section would violate unitarity since the full hadronic 
cross section for short time scales is already exploited within the constituent 
quark model. Accordingly, the 'secondaries' have to start with a vanishing 
cross section at the space-time point of the vertex or more precisely starting 
from the quark-antiquark production point that we denote by t p (cf. Fig. 3 in 
[39]) which is a fraction of 1 fm/c. During the early times, when the system is 
very dense, the 'secondaries' thus have small or vanishing cross section up to 
their hadron formation time tf (1). The latter becomes large for high p± mo- 
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menta such that the attenuation of 'secondaries' is rather low when employing 
high p±_ cuts. 
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Fig. 6. The ratio Raa (5) of charged hadrons for different centralities for Au+Au 
at y/s = 62.4 GeV. The experimental data are from Ref. [46]. The hatched bands 
correspond to our calculations with a Cronin parameter a in (4) ranging from 0.25 
to 0.4 as in Ref. [35]. 



5 Au+Au collisions at ^/s = 62.4 GeV 



The system Au+Au at y/s = 62.4 GeV has recently been investigated ex- 
perimentally at RHIC and first results on hadron attenuation have become 
available at moderate p± from the PHOBOS Collaboration [46]. We thus start 
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with single hadron attenuation in the next subsection following exactly the 
studies in Ref. [35] that have been performed at y/s = 200 GeV. Since all 
details have been presented in the former publication we step on with the ac- 
tual results. However, it is worth to point out that the 'Cronin enhancement 
factors' are higher at y/s = 62.4 GeV than at y/s = 200 GeV due to steeper 
hadron spectra in p±_ (cf. also Ref. [30]). 

5. 1 High p± hadron suppression 

The centrality dependence of the ratio Raa for charged hadrons is shown 
in Fig. 6 for 0%---6%, 6 % • • • 15 %,15 % • • • 25 %, 25%---35%, 35%---45% 
and 45 % • • • 50 % centrality of Au+Au collisions at y/s = 62 GeV. Again the 
hatched bands (color: light blue) correspond to our calculations with a Cronin 
parameter a in (4) ranging from 0.25 to 0.4 while the data stem from Ref. [46]. 
Note, that the uncertainty in the Cronin enhancement (width of the hatched 
band) decreases slightly for more peripheral reactions, which is a direct con- 
sequence of the lower number of hard NN collisions in (4). 

We emphasize, that the Cronin enhancement is most visible at momenta up to 
5 GeV/c but is practically negligible for p± > 6-8 GeV/c. As pointed out be- 
fore in Ref. [35], the suppression seen in the calculation for larger transverse 
momentum hadrons is due to the interactions of the leading (pre-)hadrons 
with target /projectile nucleons and the bulk of low momentum hadrons. Ac- 
cordingly, the ratio Raa drops below 1 in Fig. 6 in the high momentum regime 
reaching an attenuation which is about the same as that in the calculations for 
y/s = 200 GeV and p± > 5 GeV/c. Our numerical results in Fig. 6 are roughly 
compatible with the data from the PHOBOS Collaboration [46] but measure- 
ments at higher transverse momenta will be necessary to prove/disprove the 
dropping of the ratio Raa with momentum p± (cf. also Ref. [63]). 

5.2 Jet correlations from pp and Au+Au collisions 

Whereas first experimental information on the single hadron attenuation is 
available (cf. previous subsection) the quenching of 'far-side' jets at y/s = 62 A 
GeV is still unknown from the experimental side. In this respect we provide 
predictions on the basis of our HSD calculations for the angular correlations 
in pp and Au+Au collisions. 

In Fig. 7 we show the angular correlations of 'near-side' and 'far-side' charged 
hadrons by gating on high p± particles in the interval (p^™ 9 = 3 ... 6 GeVVc) 
and accumulating further charged hadrons in the interval p± = 2 GeV . . . p™ g 
< 1.0). In addition to the angular correlations from central Au+Au colli- 
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Fig. 7. 'Near-side' and 'far-side' jet-like correlations from HSD for p+p 
(dashed lines) and central Au+Au collisions (solid lines) at midrapidity for 
p T ™9 = 3 ... 6 GeV/c with p ± = 2 GeV/c . . . p]_ ri9 at yfa = 62.4 GeV. 



sions (solid lines) we have added the results from pp reactions (dashed lines) as 
resulting from the PYTHIA calculations. These PYTHIA results serve again 
as reference lines for the angular correlation without any final state interaction. 

It is seen from the pp as well as Au+Au angular correlations that the 'far- 
side' peaks are more pronounced (within the cuts taken) than the 'near-side' 
peaks. This effect can be traced back to the low probability of observing a 
further hadron close to the 'near-side' peak when gating on a high momentum 
trigger particle; most of the associated hadrons can be found in the 'far-side' 
peak. This comes about because the trigger particle is not counted in the 
'near-side' correlation whereas the 'far-side' correlation includes the leading 
particle of the partner jet. This explanation is quantified in Table 2 where the 
probability for having 0. . . 3 further charged particles in the 'near-side' and 
'far-side' region is given for pp (upper part) and central Au+Au collisions at 
y/s = 62.4 GeV. 

The relative suppression of the 'near-side' peak is very moderate whereas the 
relative suppression of the 'far-side' correlation in central Au+Au relative to 
pp reactions is ~ 60 % as for yfs = 200 GeV. This is basically due to the fact 
that the cross section of pre-hadrons with ordinary hadrons or quark-diquark 
fragments of the initial hadrons is almost independent on energy and the string 
density only moderately lower at \fs = 62.4 GeV. On the other hand, ordinary 
hadrons are allowed to form slightly earlier because the energy density (in the 
local rest frame) drops faster below 1 GeV/fm 3 . 
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Table 2 

The conditional probability to find in addition to the trigger particle 0,1,. .. 3 parti- 
cles at the far side and 0,1,. .. 3 particles at the near side in p+p and central Au+Au 
collisions at ^fs = 62 GeV. The trigger conditions are (pl ri 9 = 3...6GeV/c) and 
p± = 2 GeV . . .p T ^ 19 (\y\ < 1.0) for the associated hadrons. 

6 Summary 

Summarizing, we point out that (pre-) hadronic final state interactions con- 
tribute significantly to the high suppression effects observed in Au+Au 
collisions at RHIC. This finding is important, since the same dynamics also 
describe the hadron formation and attenuation in deep-inelastic lepton scat- 
tering off nuclei at HERMES [39,40] appreciably well. In particular, it has 
been demonstrated that the centrality dependence of the modification factor 
-Raa (5) in Au+Au collisions at y/s = 200 GeV is well described for peripheral 
and mid-central collisions on the basis of leading pre-hadron interactions [35] 
whereas the attenuation in central Au+Au collisions is noticeably underesti- 
mated. A similar observation also holds for an analysis of the 'far-side' angular 
correlations, which show a substantial, but not fully complete suppression in 
central collisions as indicated by the experimental data [9]. From these results 
one should conclude that there are some additional (and earlier) interactions 
of partons in a possibly colored medium that have not been accounted for in 
our present HSD transport studies. 

We have, furthermore, explored the possibility if shorter formations times tj 
for the hadrons might lead to an improved description of the angular cor- 
relation data since smaller Tf lead to an increase of the interaction rate at 
early times and to an increase of the elliptic flow of charged hadrons [34,62]. 
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However, for formation times t/ < 0.5 fm/c also the 'near-side' peak shows a 
reduction by the final state interactions. This reduction of the 'near-side' peak 
is still very modest for tj = 0.5 fm/c but reaches a level of ~ 50% for tj = 0.2 
fm/c which can be ruled out by the present data. In short: a sizeable reduc- 
tion of the formation time - though leading to a more substantial suppression 
of the 'far-side' peak and to an increase of the elliptic flow v-i of hadrons - is 
incompatible with the experimental observations for central Au+Au collisions. 

We close in pointing out that further experimental studies on the suppression 
of high momentum hadrons from d+Au and Au+Au collisions down to \fs 
= 20 GeV will be necessary to clearly separate initial state Cronin effects 
from final state attenuation and to disentangle the role of partons in a colored 
partonic medium from those of interacting pre-hadrons in an approximately 
color-neutral hot and dense fireball. A major problem, however, is that the jets 
escaping the medium stem from the corona of the fireball [30,31,60,64] and 
thus do not carry much information from the dense medium. Nevertheless, 
'far-side' jet suppression appears to be a promising observable since the pre- 
hadronic (color neutral) interactions employed so far lead to a suppression that 
is dominantly characterized by a linear length scale (cf. the comparison in Fig. 
4). Thus when gating on peripheral reactions with A part from 50 to 150 one 
might exploit the geometry of the fireball - which is of approximate almond 
shape - in more detail. To this aim narrow bins in A part should be performed 
and the di-hadron correlations be investigated as a function of the azimuthal 
angle of the trigger hadron relative to the reaction plane. Furthermore, also 
several cuts in the minimum transverse momentum of the associated particles 
should be taken. It is clear that very high statistics for jet correlations will 
be needed to obtain a solid information. Some experimental work along this 
direction is presently carried out. 

The authors like to acknowledge stimulating discussions with C. Greiner, T. 
Falter, B. Kopeliovich and I. Vitev throughout this study. 
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